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Abstract 

For a precise determination of the m„ and nid quark masses, it is useful to combine isospin symmetric results of lattice 
or sum rules QCD techniques with some isospin breaking study performed in chiral perturbation theory (ChPT). The 
most promising process for the later is 77 — > 3;7r decay. However, this process is affected by large chiral corrections 
and there are observed discrepancies between the values of Dalitz plot parameters stemming from the standard ChPT 
computation of its amplitude and those experimentally measured. We describe here the method based on analytic 
dispersive representation which attempts to obtain the information on the masses by taking this discrepancy into 
account independently of its exact origin. Together with the results of this analysis we present a review of some other 
constraints on the light quark masses and conclude with values of the masses compatible with them. 
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1. Introduction 

Quark masses are fundamental free parameters of the 
Standard model. Because of the QCD phenomenon of 
color confinement, quarks are bounded inside hadrons 
and their masses cannot be directly measured. In- 
stead, for their determination one uses some physical 
observable depending on the masses in a theoretical ap- 
proach and compares such theoretical prediction with 
the corresponding experimental value. Such definition 
of the masses thus depends on the theoretical approach 
used. In the following we deal with the current quark 
masses, which are the masses appearing in the QCD La- 
grangian. At low-energy region, where the color con- 
finement reigns, the non-perturbative methods of QCD 
are required. Those relevant for the determination of 
the light quark masses {niu, nid, iris) are sum rules (SR), 
lattice (LQCD) and effective field theories (EFT). 

QCD sum rules are based on dispersion relations 
stemming from analytic properties of some observable 
(such as differential decay rate of t). They connect to- 
gether the hadron world measured in experiments with 
perturbative QCD in terms of operator product expan- 
sion. For the recent review on the SR results on light 
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quark masses we refer to fP 'S^. Note that a computation 
of electromagnetic (EM) corrections within SR methods 
would be very difficult since the inclusion of long-range 
interactions changes significantly analytic properties of 
the amplitudes, which makes it impossible to write dis- 
persion relations in the regular form. 

Numerical simulations on lattice have achieved in re- 
cent years a considerable progress — the current day 
simulations are performed with 2+1 dynamical quarks, 
moreover, with near its physical value, while in order 
to reach the physical point in m„^ chiral extrapolation 
is predominantly still necessary. However, again the in- 
clusion of EM interactions on the lattice is very involved 
for their long-range character. Nevertheless, even in this 
aspect, there is some progress as quenched QED simu- 
lations are already performed (e.g. 13J). In recent years 
there appeared two attempts 10] |5l to review and aver- 
age existing lattice calculations of various quantities, in- 
cluding nig, where more details can be found. 

Both of these methods independently determined the 
values of and of isospin averaged m = "'"2"''' with 
a reasonable precision with more-or-less compatible re- 
sults. But since the isospin breaking effects on the ob- 
servables studied by both these methods generated by 
EM interactions are of the same order as those stem- 
ming from m„ - nij difference, for determination of the 
individual m„, nid they both need additional inpuj^ 

' Even for a precise determination of m EM interactions cannot be 
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A method that can provide such input is chiral per- 
turbation theory (ChPT) |6l, playing a prominent role 
among those EFT. Note that for the determination of 
the quark masses ChPT alone is insufficient. As in all 
physical results niq occurs multiplied by scalar quark 
condensate Bq and thus rescalling both of them does 
not change the physics. Moreover, starting at next- 
to-leading order (NLO) in chiral counting, there ex- 
ists a transformation of the masses [TJ, which together 
with the corresponding change in low-energy constants 
L(„Lj,Lii and some of NNLO constants C,- leave pseu- 
doscalar masses, scattering amplitudes and matrix ele- 
ments containing vector or axial vector currents invari- 
ant. Thus, using only experimental measurements we 
cannot fix this so-called Kaplan-Manohar (KM) ambi- 
guity. Consequently, ChPT can determine only quark 
mass ratios and needs some external theoretical input in 
order to fix the physical definition of the masses. For 
a review of ChPT NLO determinations of Ught quark 
masses, we refer to |8l. In addition, in |'9l the employ- 
ment of NNLO meson mass formulae is included. 

For a more precise determination of light quark 
masses it is therefore useful to combine isospin symmet- 
ric results of LQCD and SR with some isospin-breaking 
study performed in ChPT. A very suitable process for 
such study is ?; — > 3?: decay, which is possible only 
in the isospin breaking world. Moreover, EM contri- 
butions to this decay are very small and thus its ampli- 
tude is to a good approximation directly proportional to 
nid - m„. (For a more detailed introduction to this pro- 
cess see 110|.) We pull out the following normalization 
factor out of its amplitude 

v3 fii — III 

m(s,t,u) = — M(s,f,M), with R = ^ . (1) 

AR m,/ - m„ 

By computation of so defined M{s, t, u) in ChPT and 
comparing it to measured decay rate of this decay F, one 
can determine the parameter R. After pulling out this 
isospin breaking parameter we can perform isospin limit 
in the rest of the amplitude M{s, t, u), which is a very 
good (but once the experimental data on 77 — > 3;7r de- 
cays are available with a good precision, unnecessary) 
approximation. In this approximation the amplitudes 
of the neutral decay rj — > 3;r" and of the charged one 
77 -K^Ti^n^ are related together. In the following, we 
therefore work in such first order in the isospin breaking 
violation (or in other words in the limit m„± - m„i>) and 
deal mainly with the charged 77 decay. 
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b 
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0.025 + 0.160 


8 


~ 





a 


-0.030 + 0.005 


0.013 + 0.016 



Table 1 : Comparison of the experimentally measured values of Dalitz 
plot parameters with the values stemming from NNLO ChPT 



2. ij^Sn in ChPT 

Quite recently, an inclusion of the two-loop correc- 
tions to this decay amplitude was performed in [ 11 j . 

However, from the computed three successive orders 
it is obvious that chiral corrections in this process are 
large — we illustrate this on listing the values of R stem- 
ming from them and experimentally measured F |[T] : 

Rlo = 19.1, RuLO = 31.8, /?NNLO = 41.3. (2) 

Moreover, although the NNLO computation ifTTl 
leads to a reasonable result for the amplitude and for the 
value of R, there are still discrepancies between the ex- 
perimentally measured values of Dalitz parameters de- 
scribing their energy dependencies and the values pre- 
dicted from fTT^. The charged amplitude is usually pa- 
rameterized in the form (normalized to one at the center 
of Dalitz plot s = t - u = s - ^(m^ + 3m„)). 



|M,(0,0)|2 



I +ay + by^ + dx^ + fy^ + g^P'y + ... (3) 



V3(M-f) 



with variables X = 2^X^-_^'„^) , y = 2n,„(m.,-3,n,) 

The parametrization of the neutral decay reads 



\Mo(x,yf 
|Mo(0,0)|2 



= 1 + 2az + 2/3y{3z - 4/) + . 



(4) 



neglected. However, since they are here less important, at the current 
level of precision of m , it is enough to take an estimate of them. 



where in addition to the variable y, there appears dis- 
tance from the center of the Dalitz plot z = x^ +y^. 

In Table [T] we list the values of Dalitz parameters 
stemming from NNLO ChPT IfTTl and their experimen- 
tal determination from KLOE fT^l as the only measure- 
ment giving all the quoted parameters with reasonable 
precision. Note the discrepancy of the central values of 
parameters b and a. 

A third complication in the computation as well as a 
possible explanation of the discussed discrepancies in 
Dalitz parameters is the poor knowledge of the O(p^) 
low-energy constants C,. Two-loop amplitudes depend 
on subsets of 102 C,s, whose determination is needed 
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before any reliable prediction, but nowadays many of 
the C,s are just estimated (predominantly from reso- 
nance saturation). 

These complications impose the following questions: 
What is the origin of the Dalitz parameters discrepancy? 
How do this discrepancy together with the slow chiral 
convergence and the poor knowledge of the C,s influ- 
ence the determination of the isospin breaking ratio Rl 

There exist various alternative approaches taking dif- 
ferent assumptions than standard ChPT IfTSUTTl . which 
investigate further possible explanations of the discrep- 
ancy, namely higher order final state rescatterings, an in- 
fluence of a slow convergence of nn scattering or even of 
the 77 — > 37r amplitude itself. Unfortunately, the current 
situation does not enable to draw strong conclusions on 
these questions since we only have a set of 5 numbers 
from just one experiment for the charged decay and just 
one number (although independently confirmed by var- 
ious experiments) for the neutral decay. Moreover, with 
the exception of the resummed ChPT study |TT| (which 
nevertheless gives only qualitative answers) these alter- 
native approaches do not fix normalization of the am- 
plitudes (there do not appear nig explicitly). Thus, in 
order to provide any information on the quark masses, 
they unavoidably need to be matched to ChPT, which 
brings into such determination additional assumptions 
about the region where both approaches are compatible 
and give physical results. (Moreover, at the time being, 
they deal with the problem of the poor knowledge of C, s 
by using NLO ChPT results only.) 

In IfTOl we have also touched the question of possi- 
ble explanations of the Dalitz plot discrepancy but con- 
centrated more on the second question and presented a 
method that tries to extract as much information from 
the experiment and the two-loop ChPT result taking into 
account these three complications no matter what is the 
exact explanation of the Dalitz parameters discrepancy. 

For this purpose we have used our analytic dispersive 
parametrization, which is based on basic assumptions 
of QFT together with some hierarchy of various contri- 
butions to the amplitude (inspired by a very basic chiral 
counting and/or numerical studies). It takes into account 
two final-state rescatterings and can include full isospin 
violation m„± + m„o . In the first order of isospin break- 
ing the charged amplitude is parametrized in the form 

M^(s, t, u) - P(s, t, u) + U(s, t, u), (5) 
where P(s, t, u) is the polynomial part of the amplitude 

P = A,,mI + BAs -s) + CAs - sf + EAs - sf 
+ DA{t - sf + (u- sf] + FA(t - sf + (u- sf] (6) 





KLOE 1 12 1 ChPT 111] NREFT 1 16l 


reli 


0.02 ± 0. 12 -0.03 ± 0.72 0.35 ± 0. 13 


rel2 


0.12 ±0.21 -0.13 ± 1.4 0.44 ±0.20 


10' /3 


? -2 ± 25 -4.2 ± 0.7 



Table 2: Values of C; -independent combinations (CIR) of Dalitz pa- 
rameters corresponding to various determinations. 



and unitarity part U contains the parameters Ay,B^, 
Cx,Djc together with subthreshold parameters a„,f3„, 
Ai, A2 describing nn scattering (fixed from 1 18 1). 

Additional advantage of parametrization (j5]l of the 
amplitude is that by setting its parameters to their partic- 
ular values, it can reproduce the NNLO ChPT amplitude 
exactly (on the region below the nr] thresholds). 

We have presented two analyses using different as- 
sumptions on the physical amplitude trying to use the 
information about it stemming from NNLO ChPT to- 
gether with the one from KLOE. As we have already 
noted, since the original data sets of KLOE are inac- 
cessible, we had to construct a distribution from the val- 
ues of five Dalitz parameters published by them (includ- 
ing their errors) and all our current analyses depend on 
the assumption that such distribution describes well the 
genuine physical amplitude, thereby relying on KLOE 
error determination. 

3. First analysis: Correcting ChPT rj -* 2>7i result 

Our first analysis is inspired by the following obser- 
vation. At the two-loop level, one can find some re- 
lations between Dalitz parameters that are independent 
on the values of the C,s. They read 

reli:(4(fe-H£/)-fl2-16Q')|^ = 

rel2: {a^ - Aab + Aad + 8/ - 8g)|^ = (CIR) 

rel3:jS|c = 0. 

If we compare the values of these combinations coming 
from KLOE with those from NNLO ChPT computation, 
we find a good correspondence between their central 
values, as is obvious from Table [2] (For comparison we 
have also included there the values corresponding to the 
alternative approach of 1 16| (NREFT), where a possible 
explanation for the a discrepancy is proposed.) 

It motivates us to assume in this analysis that all the 
discrepancy can be included into a small real polyno- 
mial contribution (i.e. into a possible change of the val- 
ues of the C,s). 

Before we proceed to such analysis let us make a few 
comments. The error bars quoted in Table|2]are overes- 
timated since they were obtained just as combinations of 
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/>(4) 


0.46(1) 


1.95(10) 


-0.6(2) 


1.04(2) 






/>(6) 


0.58(1) 


2.4(2) 


0.3(34) 


1.6(24) 


5(150) 


-4(84) 


p(6) 
^ corr. 


0.575(6) 


1.99(4) 


-6.8(3) 


0.94(3) 


-31(3) 


20(1) 



Table 3: The values of parameters of Jsj corresponding to the "cor- 
rected amplitude" of the first analysis. 



the errors of the individual parameters — this is the case 
especially for the ChPT values, where the parameters 
were affected by the large C, uncertainties, whose effect 
on these combinations is reduced. Hence, it would be 
desirable to remeasure these combinations (and reeval- 
uate them also in the ChPT analysis) in order to test this 
class of explanations. A case of a special importance 
is the neutral parameter /3, which is thus determined in 
ChPT independently on the C,s. 

Once this correspondence is confirmed, from the val- 
ues of the individual Dalitz parameters [or more suitably 
of those of parametrization (j5])] corresponding to the ex- 
perimental amplitude one could construct sum rules for 
the C,s and examine their compatibility with the values 
coming from other observables. 

Now let us return to our analysis. We assume that 
by an addition of a small real polynomial AP to the 
NNLO ChPT amplitude one would reproduce the phys- 
ical data. By fitting such AP from KLOE we obtain 
such "corrected amplitude". We have found IITOll that 
the corrections are indeed small (on the physical region) 
and correspond to the change of parameters in polyno- 
mial part (j6]), which is shown in Table [3] Even without 
changing its unitarity part, Mcon. reproduces very well 
the KLOE-like distribution and the corresponding value 
of R is shifted to 

R = 37.7 + 2.9 [ChPT-t-Disp.-HKLOE]. (7) 

The quoted error is estimated conservatively taking into 
account the slow convergence of chiral expansion in the 
first three orders — assuming that the error of the cor- 
rected amplitude is equal to one half of the difference 
between this amplitude and the NLO one (similarly as 
it was done with the NNLO ampUtude in ifTTll '). 

4. Second analysis: Direct fit to // -» 3;r data 

A different analysis not imposing anything about the 
origin of the Dalitz parameters discrepancy takes advan- 
tage of the fact that the construction of parametrization 
(|5]l employs just general properties of the amplitude. We 
can thus assume that the physical amplitude (for now the 
KLOE-like distribution) should be reproduced by this 
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Figure 1: Comparison of the amplitude obtained by our second 
analysis (solid lines) with the order-by-order parametrization of the 
NNLO ChPT amplitude (dotted) and its "resummed" parametrization 
(dashed) on the cut t = u. (More details in the main text.) The vertical 
lines indicate the physical region. 

parametrization, which gives us a clear prescription for 
its fitting. In order to fix the normalization we suppose 
that ChPT determination of the amplitude is reliable at 
least in the regiorj^ specified below. 

The usually employed fact when searching for such 
matching region is the accidental coincidence at the chi- 
ral NLO level of the following points (on the cut s - u): 

a) the point, where the amplitude is of order 0(m^) 
(Adler zero), 

b) the point, where Re M = 0, 

c) the point, where the corrections of the computed 
order to the slope are small. 

However, this coincidence proves not to be the case at 
NNLO. Therefore, having the NNLO results at hand 
there is no reason why the amplitude at e.g. the point 
fulfilling condition b) should have a faster chiral con- 
vergence than at any other point (in this particular case 
at NNLO it even seems to be the opposite). 

Instead, taking parametrization (|5]l and the first three 
chiral orders of the amplitude that are at disposal, we 
have found a more suitable prescriptiorj^ [10 1: We 
should fit on the t - u cut and match only imaginary 
parts of the amplitudes under the physical thresholds 
(i.e. the dependence on the C,s is reduced). Moreover, 
one can further reduce the error of this procedure by fix- 
ing the normalization so that the resulting parametriza- 
tion in this region interpolates between the order-by- 
order and the "resummed" fits (|5]l of the chiral am- 
plitude (they differ in the unitarity part, in the choice. 



^Note that since we have analytic results, for fixing of the normal- 
ization one matching point would be sufficient. 

■'However, it could happen that in future e.g. some NNNLO com- 
putation reveals that even this prescription is not the best one. 
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whether one respects chiral orders of the parameters ap- 
pearing there or whether one reorders them taking into 
account also some higher unitarity contributions). 

In Figure [T] we plot the result of this matching in 
comparison to the order-by-order and the "resummed" 
parametrization of the original NNLO ChPT amplitude. 

By comparing measured T 1 1 1 with integration of the 
resulting amplitude over phase space, we obtain 



R = 37.8 + 3.3 [Disp.-HKLOE]. 



(8) 



The possible sources of errors entering this analysis in- 
clude: the uncertainties of the experimental data to- 
gether with the (in)accuracy of their description by our 
parametrization; the errors induced by analytic contin- 
uation of the parametrization to the region where we 
match to ChPT; and finally the error of chiral expansion 
of the amplitude (its 0{p^) order) which is used for the 
matching in this region. Taking into account properties 
of parametrization (|5]l and the above described construc- 
tion of the matching procedure, there is no surprise that 
the first type of errors prevails over other two even if we 
estimate them conservatively. 

Further details on both analyses can be found in ifTOll . 

5. Results 

Since the dominant errors in both the performed anal- 
yses are of diff'erent origin (and the values are compati- 
ble), we can combine them into 



R = 37.7 ± 2.2. 



(9) 



However, as was stated above, both analyses rely on the 
assumption that the genuine physical amplitude is de- 
scribed by the distribution constructed from the values 
of Dalitz parameters given by KLOE (and on the fact 
that one can use the NNLO ChPT result at least as the 
input for the matching). Without further measurements, 
the error bar connected with this assumption is difficult 
to quantify. We estimate it from the second analysis and 
use in the following also more conservative value 



R = 37.7 ± 3.3. 



(10) 



In Figure [2] we review various constraints on the 
light quark masses coming from various recent analy- 
ses. In order to keep the plot uncluttered, we include 
only those which are independent and when diff'erent 
updates of some analysis exist, we list only the most re- 
cent one. For older results we refer to |8|. The axes 
of the presented plot correspond to (isospin symmetric) 
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Figure 2: Various constraints on quark mass ratios (description given 
in the main text) 



ratio r = ^ and to quantitative measure of isospin vio- 
lation which appear naturally in various analysesj^ 

Our result (j9]l is denoted by Ri whereas the more con- 
servative one by R2. The results of the alternative analy- 
sis ofT]—>3n from lfT4l[T5l are depicted as Qi (with the 
error bar estimated by fS\) and the more recent resulj^ 
with Q2. This numeric approach uses dispersive rela- 
tions of Omnes type, thereby attempting to include two- 
pion rescatterings to all orders. The necessary matching 
to ChPT is performed at the NLO Adler zero (see above) 
to the NLO amplitude. For r > 26 it leads (similarly as 
ChPT NLO computation) to m„ lighter than our result. 

Interestingly, our result is fully compatible with the 
result of 15] using ChPT NNLO expressions for me- 
son masses that take into account only the EM contri- 
butions due to Dashen |fT9] (central value and error bars 
denoted by Mi). The result of |9| that includes the vio- 
lation of Dashen limit computed in |20| xd - 1-84 (de- 
picted without error bars as M2) moves m„ to noticeable 
smaller values. If one used only NLO expressions for 



^The different plot used in previous summaries (as e.g. fE\) is con- 

2 - 2 



nected witli anotlier isospin breaking parameter Q 



— whicli 



liave at NLO in CliPT the advantages that it can be expressed in terms 
of QCD pseudoscalar masses only and is reasonably stable with re- 
spect to KM transformation. However, at the current level of precision 
one includes also chiral two-loop eifects. with which both of these ad- 
vantages of Q are lost since the relation between Q and the meson 
masses gains noticeable r-dependent corrections at NNLO (cf. |9|). 

Also in our analyses, instead of M we could have used Mq defined 
by pulling out . However, since we need to match to NNLO ChPT 
nil , the normalization with R is more natural. Note that KM ambi- 
guity is assumed to be fixed in 1 1 1 J by the values of LLCs used there 
(e.g. the value of = stemming from large considerations). 
'Talk by E. Passemar at EPS HEP 201 1, Grenoble, July 201 1. 
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the masses, for the value of r < 28 our result would lie 
in between the determination respecting Dashen (M3) 
and the one with xq = 1.84 (slightly above the upper 
bound of 22)- For r < 27 this NLO determination vio- 
lating Dashen limit would be more compatible with our 
result than the one respecting it. 

We proceed now to the recent determinations of r. 
Progress report on sum rules [T| quotes the value de- 
noted by r4. The method of [2 1 J inspired by large Nc 
considerations enables to determine r from the ratio of 
r( J'Zyy') ■ Ths original number is indicated by arrow 
'■5. In II22I such analysis was repeated with more re- 
cent numbers, giving the only recent determination of r 
incompatible with the lattice averages. There was also 
used the result of Ii23il for Dashen violation x^, = 1.5 and 
from niK* - m^o obtained a different value of Q. These 
two results of (22] together are depicted with label KN. 

The most advanced determinations of r are performed 
using lattice. Nowadays there exist two averaging at- 
tempts on lattice results. The average of Laiho et al. ^ 
gives ri = 27.55(14) and isospin quark masse^ 

n/, = 93.6(1. l)GeV, m'' = 3.419(47) GeV. (11) 

FLAG group [4| have similar averages but in order to be 
more conservative, they quote estimates ry, = 27.4(4), 



= 94(3) GeV, mf = 3.43(11) GeV. 



(12) 



We use FLAG numbers as our conservative choice. 
However, since their estimate are in some aspects very 
strict (e.g. not including BMW |24| yet), in order to 
show the precision of the current determinations we also 
use the average of |5 1 with error bars extended so that 
they include all the central values of the individual re- 
cent lattice determinations, r2 - 27.55(25), 

= 93.6(2.2) GeV, m^' = 3.42(9) GeV. (13) 

In Figure |2] we have also included recent isospin 
breaking study on lattice O denoted by B (note that it 
is included in average [31 but not in f4l). For complete- 
ness, let us note that baryonic determinations as 
give smaller quark ratio < 0.4. 

Taking our result 



mate ( [T2| ), we obtain 

m.. 

0.48(3) [0.48(2)]^ 



10 1 together with the lattice esti- 



2 = 23(1) [23.2(7)], 



MS,//=2GeV 



= 2.23(13)MeV [2.21(9) MeV] , 
m^^'^^^'^"'' =4.63(16)MeV [4.62(12) MeV] . 



The more precise values given in the square brackets 
correspond to the situation if our assumptions on the 
77 — » 3;r data and the estimates on the isospin quark 
masses ( [T3| l were fulfilled. 
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